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1 .  Introduction 


The  study  of  time-dependent  oscillator  systems  has  received  a  great  deal 

of  attention  both  in  classical  and  quantum  mechanical  studies.  In  the  usual 

quantum-mechanical  treatment  of  these  systems,  a  time  -  independent  Hamiltonian 

is  assumed,  which  must  be  obtained  through  the  Legendre  transformation  of  the 

Lagrangian  function  which  describes  the  equation  of  motion  for  the  system. 

In  this  context,  it  is  unlikely  that  all  physical  oscillator  systems  have 

conserved  Hamiltonians .  The  Hamiltonian  of  a  damped  driven  harmonic 

oscillator  is  a  typical  one  with  explicit  time-dependence.  Even  though  there 

have  been  extensive  worksv  '  on  this  system,  this  apparentlv  simple  system  has 

not  been  exactly  solved  quantum  mechanically,  in  contrast  to  the  classical 

case.  Furthermore,  the  physical  interpretations  of  the  Hamiltonian  itself  and 

of  the  time -dependent  Schrodinger  equation  obtained  by  replacing  the  canonical 

momentum  with  ?  —  in  the  Kaldirola-Kanai -  type  Hamiltonian^  ’  differ  among 

(4) 

researchers.  Some  claim  that  this  Hamiltonian  (accordingly  the  Schrodinger 
equation)  describes  a  dissipative  quantum  system,  while  others object  to 
such  an  interpretation.  From  the  viewpoint  of  group  theory,  Cervero  et  al^ 
interpreted  it  to  represent  a  variable  frequency  oscillator  with  minimum 
uncertainty . 

The  Feynman  path  integral  formalism*^  provides  a  general  approach  to 

/  o  \ 

quantum  systems.  Montroll'  transformed  the  original  Feynman  path  integral 

(9) 

into  a  Gaussian  integral,  and  his  method  was  applied  by  Cheng  to  evaluate 
the  propagator  for  a  forced  harmonic  system.  Once  the  propagator  suitable  to 
a  given  system  is  obtained,  the  wavefunction  of  the  system  is  easily 
calculated.  However,  even  if  the  Schrodinger  equation  is  solved  through  the 
path  integral  or  some  other  method,  there  remains  the  problem  of  second- 
quantization.  In  general,  the  Hamiltonian  for  the  damped  driven  harmonic 


oscillator  (DDHO)  is  not  conserved  and  not  identical  to  the  total  energy  of 
the  system,  giving  rise  to  much  difficutly  for  developing  a  second- 
quantization  formalism.  However,  second-quantization,  if  possible,  is 
important  in  connection  with  the  problem  of  finding  invariants  of  a  non- 
conserved  system,  and  for  obtaining  the  coherent  states of  DDHO. 

On  the  other  hand,  the  optical  properties  of  molecules  are  greatly 

altered  when  they  are  adsorbed  cn  or  near  a  solid  surface.  During  the  past 

decade,  many  theoretical  and  experimental  methods  have  been  developed  to 

investigate  surface  spectroscopy,  surface  -  enhanced  Ramann  scattering 

(SERS),^^  resonance  fluorescence ,  '  photodissociation, etc.  There 

exist,  however,  several  criticisms  and  controversies  regarding  the 

(14) 

electromagnetic  theory.  In  order  to  better  understand  the  surface 

phenomena  mentioned  above,  one  must  study  at  the  quantum  mechanical  level  the 
dynamical  behavior  and  optical  response  of  the  molecules  adsorbed  on  a  solid 
(metal  or  dielectric)  surface. 

In  this  work,  we  present  a  study  of  the  quantum  mechanical  properties  of 
a  molecule  adsorbed  on  or  near  a  dielectric  plane  surface.  We  represent  the 
molecule  as  a  polarizable  point  particle  subject  to  a  time-dependent  external 
electric  field.  We  do  not  consider  the  nuclear  vibrations  in  this  work.  In 
Sec.  2  the  classical  equation  of  motion  of  the  induced  molecular  dipole  is 
given,  and  this  equation  is  transformed  into  that  of  a  damped  harmonic 
oscillator  driven  by  an  external  force.  In  Sec.  3  the  propagator  and 
wavefunction  are  obtained  through  the  Feynman  path  integral  method.  Energy 
expectation  values  and  dipole  moment  matrix  elements  are  expressed  in  Sec.  4. 
See.  *  is  devoted  to  second  quantization  and  calculation  of  the  transition 
probability  between  the  harmonic  oscillator  states  and  damped  driven  harmonic 


oscillator  states.  In  Sec.  6  we  discuss  the  results. 


2 .  Classical  Equation  of  Motion  for  the  Dipole  Moment 

We  consider  a  molecule  adsorbed  on  a  dielectric  plane  surface  with 

complex  dielectric  constant  e  and  light  beam  incident  from  the  vacuum  side  at 

an  angle  S  with  respect  to  the  z-axis  normal  to  the  surface  (xy-plane) .  The 

normal  modes  of  the  incident  field  satisfy  Maxwell's  equations  subject  to 

appropriate  boundary  conditions.  The  wave  vector  k  of  the  field,  which  is 

assumed  to  be  independent  of  frequency,  is  decomposed  into  a  two-dimensional 

component  ki|  parallel  to  the  surface  and  a  z -component  perpendicular  to  the 

surface.  The  electric  field  vector  is  polarized  in  two  ways.  The  first  is 

the  s-polarized  wave  (s-wave)  perpendicular  to  the  incident  plane  defined  by  k 

and  z  (a  vector  along  the  z -direction) .  The  second  is  the  p-polarized  wave 

(p-wave)  in  the  incident  plane.  The  local  field  E_(r  ,t)  at  the  adsorbed 

i  m 

—¥ 

molecule  position  r^  is  expressed  as 

ZAvm  ,t)  -  E  (r  t)  +  E  (r  t)  .  (2-1) 

if  m  pm  s  m 

— ♦  — ► 

where  t  indicates  time.  E^r^.t)  represents  the  primary  field,  which  consists 

of  the  field  incident  directly  on  the  molecule  and  the  field  reflected  from 

— ¥ 

the  solid  surface  when  the  molecule  is  absent,  and  E  is  calculated  using  the 

P 

(15)  -  - 

Fresnel  formula.  E  (r  , t)  is  the  secondary  field  which  represents  the 

field  emitted  by  the  induced  molecular  dipole  and  reflected  back  to  it  from 
the  surface.  The  induced  dipole  moment  p(t)  is  the  sum  of  dipole  moments 
induced  by  the  primary  and  secondary  fields.  The  classical  equation  of  motion 
for  Ai(t)  is  given  by 


M  +  7P  + 


‘V 


wQa. 


(2-2) 
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The  docs  indicate  cime  derivatives,  u>q  is  the  oscillation  frequency  of  the 
dipole  charge,  7q  is  a  natural  damping  constant  which  measures  the  linewidth 
in  the  optical  processes ,  and  a  represent  the  effective  molecular 
polarizability  tensor  of  second  rank.  In  Eq .  (2-2)  we  neglected  the  magnetic 
effects  on  the  dipole  motion,  and  hereafter  we  set  the  permeability  of  the 
surface  dielectric  to  be  unity.  At  the  molecular  position  (0,0, d)  the  primary 
field  vector  is  calculated  using  the  Fresnel  formulas  as 


E 

P 


(d,t)  -  Eq 


exp 


-i(kd  cosB 


+  u' t ) 


(2-3) 


where  the  components  of  Eq  are  given  as 


E0x  " 

E0,  -  Et>{1  + 

EQz  -  Eo{l  +  f^(«  (u>'  )  ,  8)  j'sintf 
f ,<<(„■>,«)  - 

t (w' )  +  [ e (w' )  -  sin  9  )  2 


exp(2ikdcos0 ) 


f  2(e(u>'),8) 


2  4 

COS0  -  ;  £  (u>'  )  -  sin  8  1 

“ — *T  exp(2rkdcos0 ) 

cos 8  +  [  «  (cj'  )  -  sin  0  ] 


e (w' )  -  €^(w' )  +  i«2(w' ) 


(2-4) 


(2-5) 


(2-6) 


Here  k  -  |k|  is  the  magnitude  of  the  wavevector,  and  w'  is  the  angular 
frequency  of  the  incident  field. 


The  secondary  field  at  the  molecule  position  is  expressed  as  a  function 
of  the  induced  dipole  moment, 


Es(d,t) 


Gs[d;d;t-t'] 


dMi-t') 
dt ' 


(2-7) 


«.  (17  18) 

where  is  the  scattering  part  of  the  dyadic  Green's  function. 

Equation  (2-7)  represents  the  self -polarization  effect  of  the  induced  dipole. 

The  electric  field  emitted  by  the  dipole  at  an  earlier  time  t'  interacts  with 

the  surface  atom  (or  molecule) ,  causing  it  to  emit  lights  by  polarization. 

This  emitted  field  by  the  surface  atom  (molecule)  interacts  with  the  adsorbed 

molecule  at  a  later  time  t  >  t' ,  adding  an  extra  induced  dipole  moment,  which 

emits  radiation  that  polarizes  the  surface  atom  (molecule)  again,  etc. 

Through  these  self -polarization  processes,  the  secondary  field  yields  a  change 

/ION 

in  the  linewidth  and  oscillation  frequency  of  the  dipole.  Efrima  et  al^ 
used  a  Fourier  transformation  method  to  solve  Eq.  (2-7),  which  gives  an 
integro-dif ferential  equation,  and  applied  it  to  the  Raman  scattering  problem 
of  an  adsorbed  molecular  system.  The  Fourier  transform  of  the  dyadic  Green's 
function  Gs(d;d;w')  can  be  calculated  only  numerically,  since  it  is  given  in 
complicate  integral  form.  To  avoid  such  complexity,  we  incorporate  the 
effects  of  the  secondary  field  into  the  linewidth  and  the  oscillation 
frequency  of  the  dipole  moment.  In  this  manner,  we  obtain  a  modified 
classical  equation  of  motion  for  the  dipole  moment, 


-♦  2-*  “ ► 

M  +  7M  +  -  cjQa(a>'  )  #Ep(d,  t) 


(2-8) 


7 


As  mentioned  above,  the  self -polarization  effects  are  absorbed  within  the 
modified  damping  constant  7  and  modified  frequency  w  whose  explicit  forms 


+  lr 


2  2 
w  =  (wn  +  Aw  ) 
m  0  m 


wZ 

=  —  n«^a(w' ) • [ Im  G^ (d ; d ; w' ) *n] ^ 


(2-9) 


Aw  -  w„ 
m  0 


-  h*-|a(w'  )  •  [Re  G^  (d;  d;  ,  w' )  «hj-j 


where  G(d;d;w')  is  the  Fourier  transform  of  Gs(d;d;t-t')  and  n  -  /i/ 1 m I  - 
Let  us  now  introduce  a  coordinate  system  in  which  the  effective 
polarizability  tensor  a  is  diagonal. 


u(w') 


a  (w' )  xx  +  a  (w' )yy  +  a  (w')zz 
x  v  J  J  Z 


(2-10) 


where  x,  y  and  z  are  unit  vectors  along  the  coordinate  axes.  We  assume  that 
the  molecule  is  isotropic  and  its  dipole  moment  is  directed  toward  the 
positive  z-direction.  In  this  case  we  can  set  n  =  z,  and  thus  only  the  z- 
component  of  the  primary  field,  which  is  given  by  Eq .  (2-3)  through  (2-6),  is 
required  in  the  equation  of  motion.  Other  molecular  dipole  moment 
configurations  are  similar  to  this  case,  so  that  our  method  does  not  lose 
generality.  By  this  manner,  Eq.  (2-8)  can  be  rewritten  as 


•  • 


(2-11) 


M  +  7P  +  w2^  -  f3(€(w'),9)  e  iw 


where 


f3(e(w' ) ,6) 


=  EQa(a/  )^gS  intf^l  +  f  ^  (  £  ( w '  )  ,  5  )  |e  '  i 


ikd  cos& 


(2-12) 


Equation  (2-11)  can  be  written  in  a  more  familiar  form  if  we  take  just  the 
real  part  of  the  external  driving  force.  After  some  calculations,  we  obtain 
the  classical  equation  of  motion  for  the  dipole  moment  of  the  adsorbed 
molecule  as 


X  +  7X  +  ii)  X 
m 


f (t)/m 


0 


(2-13) 


Here  Qx  -  p,  where  Q  is  the  average  dipole  charge,  and  raQ  is  the  mass  of  the 
molecule.  The  form  of  equation  of  motion  (2-13)  is  the  same  as  that  of  the 
usual  damped  driven  harmonic  oscillator,  so  that  f  is  given  as 


f(t)  -  g(c o'  ,6)  cos(cj't -4>  ) 


(2-14-a) 


g( 


2E0m0w0  sin*  f  2  2  j 

"  - Q~g:(—  7) - {ft”''*)  + 


g ^(w',0)  -  cosd  cos (kdcostf  )  |  e  (w' )  |  2  cosS  +  R(w'  ,  0)-|e^(w'  )cosr] 

){«l(w'  )sinr7 


+  £2(w' )sinr? 
-  e  2  (w> )  cost; 


+  R (a)', 0)  cosff  sin(kdcos0 )< 


g2(w',5)  -  R(u',ff)  cosff  cos(kdcos$  )•^£2(u>,  )cos,7  -  « ^  (<w' )  sinr;|’ 


-  R(w',0)  sin(kdcos0) 


R(w^0)  +  coss|c^(w'  )cosr;  +  t  0  (u)'  )  sinr^j 


2  o  ?  f 

g^(u'ff)  -  R“(w',0)  +  I  «  (  oj  '  )  {  cos  ^  (9  +  2R(w',0)  cosH  e  ^  (w '  )  co 


sv 


+  e  -(u1  )  sinp 


} 


(2-14-b) 


<pQ  -  tan 


1|g2 ' €  (w'  )  .  ^)/g1(t  («*>' )  ,  0)  j 


T7  -  7  tan  1|e2(w' )/[  €^(o>' )  -  sin20]| 


(2-14-c) 


R(w' ,8) 


{[VW> 


2  2  )  k 

)  -  sin  8 ]  +  £ j (w< ' 


->}' 


To  avoid  the  complicated  procedure  of  evaluating  the  wavefunction  and  other 
quantities  involved,  we  have  taken  the  real  part  of  the  driving  force.  The 
classical  Lagrangian  and  Hamiltonian  function  corresponding  to  the  equation  of 
motion  (2-13)  are  given,  respectively,  as 


T  ,  [1*2  1  22  1  ,  . 
L(x,x,t)  =  j^m^x  -  JiHqU^x  +  exp  [  yt  ] 


H(x , p , t)  - 


r2—  e  2yt  +  ^m.w2x2  -  f(t)x 
2m^  2  0m 


exp(yt) 


(2-15) 


(2-16) 


where  p  is  the  canonical  momentum. 


3 .  Propagator  and  Wavefunction 

Recently  Um  et  al^^  have  treated  a  quantum  damped  driven  harmonic 


oscillator  with  the  path  integral  method, 
solve  the  Schrodinger  wave  equation 


(7) 


Following  their  formalism,  we 


iX  ^  =  riu>(x.t 

0  t 


(3-1) 


where  the  Hamiltonian  operator  fl  is  obtained  bv  making  the  replacement  p  -  iX 

n 

in  Eq .  (2-16).  Adopting  a  Gaussian  form  for  the  propagator  K, 


JX 


Kfx(t) , t , x ( 0 ) ,0 


Jt 


.7t/J 


=  Anexp|-a(t)x  e  -  b(t)xe  -  c(t) 


(3-2) 


which  satisfies  the  wave  equation 


(3-3) 


we  obtain  three  first-order  differential  equations  for  the  coefficents  in  the 

propagator.  These  differential  equations  are  readily  solved  by  the  boundary 

condition  that  K  reduces  to  the  propagator  for  a  damped  harmonic 

(21) 

ocsillator  when  there  is  no  driving  force.  Through  the  above  procedure  we 
obtain  an  explicit  form  for  K, 


7t/2 


K(x(t) , t ; Xq , 0 j  - 


r  raowe  1 

4  1 

2 

_2ni  sin(wt) 

exp  j 

A1(t)x  -  A2(t)x 

2 

x  exp  -A  (t)xn  +  A  (t)x 


(3-4) 


11 


.2  2  , 4 

where  x^  -  x(0)  and  w  =  -  7"/4j  '  represents  the  reduced  frequency,  which 
will  be  assumed  to  be  real  throughout  the  present  text.  The  coefficients  in 
Eq .  i  3-41  are  given  as 


l(w'  .  9  )  -  g(w'  ,  $  )  /b^|  (w'  ) 


b^ (w' 1 


(7 w '  )  *"  +  ■  w'*"  +  7^/4  | 


.^2  .  f,2  ..,2  ,  2//.\2lH 


-if  2  2  2  1 

-  tan  ^7«'/(<*>  -  w'  +7/4)1 


IU_UJ  /• 

Al('}  =  IT*  e7  {21  -  cot(wt) 


} 


A2(t)  -  e~rC  [{2I"  coc(wt)|"  cos(w' t-^0-<^1)  -  ~  sin(o>' t-<^o-01) 


(3-5) 


A3(t) 


JL_  Jx_ 

2iY  \2w 


+  COt  (tilt  ) 


} 


A(w' , S )w 
A4<'t')  iW  sin(ut) 


.  /  <-\  ITIaCOX  ,  r\ 

7t/2  ......  0  -yt/2 

e  cos(w' .)  +  — - : — 7— —  e 

0  1  in  sin  ut 


The  wavefunction ,  i/>  (x,t),  of  the  dipole  moment  at  time  t  is  obtained  bv 

n 


the  formula 


dxQ  K[x,  t ;  Xq  ,  0  ]  V>n(xQ,0)  ,  (3-6) 

where  i/i^Xq.O)  is  the  wavefunction  of  a  damped  harmonic  oscillator  with  the 
reduced  frequency  w.  Combining  Eqs.  (3-4)  and  (3-6),  we  obtain 


V>n(x,  t) 


r 
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From  now  on,  we  shall  use  the  notation  £(t)  =  ►,  ►  (t)  -  f^...  D(t)  -  D,  etc. 
for  simiplicity  whenever  there  is  no  ambiguity.  In  obtaining  the  wavefunction 
we  nave  used  the  identity 


exp 


cot 


r 

L 


0 

1  +  v 


.  ,  n 
i] 


(  3  .  Q ) 


We  find  the  following  useful  relations  between  the  coefficients  in  the 
exponential  terms  of  the  wavefunction, 

B1  +  B*  =  D“.  30  +  B?"  =  2D2E,  B3  +  B*  =  D2E2  ,  (3-10) 

which  are  used  in  later  calculations. 

4 .  Energy  Expectation  Values  and  Matrix  Elements  of  the  Dipole  Moment 

Using  the  Lagrangian  of  Eq .  (2-15),  the  mechanical  energy  operator  of  the 
system  is  easily  obtained  as 


E 

op 


1  2  2 
+  n  mnw  x 
2  Om 


xf(t) 


(4-1) 


We  note  that  above  energy  operator  is  not  identical  to  the  Hamiltonian  ft 

itself.  ft  does  not  represent  the  total  energy  of  the  system,  but  rather  is 

(22) 

the  generator  of  the  motion  of  an  energy-dissipative  open  system. 

Equation  (4-1)  is  the  same  as  Eq .  (6-1)  in  Ref.  20  except  for  the  xf(t)  term. 
It  is  straightforward  to  calculate  the  matrix  elements  oi  the  energy  operator: 


<E  > 


0  m 


op  m .  n 


2D^  L 


V\n+1)  (n+2)  (t-  -  2A,B,)6  „  +  72(n+l) 

z  il  ra,n +  z 


x  (DE  +  2\. A-B. -  DA. ) 5  .  + 

1  z  1  J  m,n+i 


{<n+2 


)  +  2(2n+l)A1A4B1  +  D2E2  -  A.^2  -  2D2EA3M 


,  725  <DE  -  2JjA2A4  -  DA,)^  ^  +  -  2A1A^/)6m  n  2 


2^2 , 


x  exp 


i(m-n)  cot  ^ ( f ( t ) ) 


(4-2) 


where 


A1(C) 


m  1 


A2(t)  -=  D[B2  -  2B1E] 


A3(t) 


A («'  J)hl(u>'  ) 


ni-u 
0  m 


2  cos(w't -<f>Q)  ,  A4(t)  -  D  -  B1 


(4-3) 


Evidently,  only  the  diagonal  element  <Eop>n  n  and  four  off-diagonal  elements 


<E  >  ,,  and  <E  >  ,0  have  non-zero  values.  The  matrix  elements  of  the 
op  n±l,n  op  n±2,n 


dipole  moment  are  expressed  as 


Q_ 

<M>  -  oT:  exp 

^  m,n  2D  r 


-1. 


i(m-n)  cot  (f) 


72 (n+1)  S  .+  J2n  6  ,  +  2DE  5 

m,n+l  m,n-l  m,n 


(4-4) 


This  equation  contains  the  diagonal  element,  which  in  fact  is  closely  related 
to  the  parity  problem  of  the  wavefunction ,  as  will  be  explained  later. 

However,  the  selection  rule  for  the  dipole  transition  between  different  states 


is  Am  -  ±1.  Combining  Eqs .  (3-9)  and  (4-2),  we  arrive  at  the  following 
explicit  form  for  the  dipole  matrix  elements: 


<u>  ,  -  Q 

n+1 ,  n  v 


rLn+llWe7t:' 

4 

iwt  y  .  . 

l 

3 

O 

a 

On] 

e  +  ~  sin(ut) 

=  </i>*  . 

n,n+l 


<M>n,n  "  °A(m0'n  cos(“'  (t  "  0) 


(4-5) 


<fi>  -  0  (t  -  0) 

n ,  n 


5 .  Second  Quantization  and  Transition  Probability 

To  obtain  a  second-quantization  formalism  for  a  molecular  system  adsorbed 
on  a  solid  surface,  we  define  two  functions, 


u(t) 


?(t)  +  i 
D( t)72 [ l+£2 ( t) ] 


(5-1) 


rKt)  -  2  ij<u  (  t )  B1(t) 


where  D(t)  and  £(t)  have  already  defined  in  Eq .  (3-8)  (rj(t)  is  not  the  same  as 
the  angle  rj  in  Eq .  (2-14-b)).  Using  the  relation  of  Eq .  (3-10),  we  obtain 


T(  t)  u*(t)  -  u  ( t )  t)*(  t)  =  i# 


(5-2) 


Thus  we  define  the  time-dependent  annihilation  and  creation  operators  as 


1 


a'(t)  -  ii 


{ 


u*(t)  ^p  +  i>([B, 


-2B1E] 


rj*(t)(x-E(t)) 


(5-3) 


where  B^(t)  ,  BjCt)  and  E(t)  are  defined  in  Eq .  (3-8).  Ue  note  that  B^(t)  - 
2B^(t)  E(t)  is  pure  imaginary.  We  shall  hereafter  use  the  notation  of  u(t)  - 
u  and  rj(  t)  -  tj  .  It  is  obvious  that  the  non-Hermitian  operators  a^(t)  and  a(t) 
satisfy  the  commutation  relation 


[ a( t) ,  a1 (t) ]  -  1 


(5-4) 


By  the  definition  of  Eq .  (5-3),  we  get 


x-u*  a(t)  +  u  a ' ( t )  +  E 


p  -  r>*  a(t)  +  tj  a^(t)  -  i>t [ B0  -  2B^I 


Obviously  x  and  p  satisfy  the  relation  [x,p]  -  ij<.  Now  we  set 


|m(t)>  -  y==  (a!(t))m|0(t)> 


a( t) j  0 ( t )>  -  0 


(5-5) 


(5-6) 


(5-7) 


Combining  Eqs .  (5-3)  and  Eq.  (5-7),  we  obtain  the  time -dependent  ground-state 
wave  function 


0o(x,t)  =  <xjO(t)> 
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rmowe  -7t)  -  ii 

9 

wW1(t)  i(t)  +  i 

exp 

-B1(t)x"  +  Bj ( t)x  -  B3(t) 

(5-8) 


We  note  chat  0Q(x,t)  reduces  to  V>q(x)  at  t  -  0 ,  which  is  the  ground-state 
wavefunction  of  a  simple  harmonic  oscillator.  Combining  Eqs .  (4-1),  (5-4)  and 
(5-5),  we  find  the  second-quantized  form  of  the  energy  Hamiltonian  operator, 
respectively,  to  be 


op 


rr;2  -  2  y  t  1 


2  2] 
u 

0  m  j 


a^(t)a' (t)  + 


*_  -  2yt  1  2  2 

—  e  +  rm„w  u* 


2m 


2111  /-.w  V. 

0  ra 


a(t)a(t) 


n  -  2yt  2.  1 2]  \ 

*  *  e  +  mno)  u  (a1  ( 

m„  0  m1  '  j 


t)a(t)+-)  + 


^a(B2-2B1E)e"27t  -  uf 


m 


0 


+  m_cj  uE 
0  m 


a!(t)  + 

-  ^-(B  -2B  E)2e‘27t  +  innw2E2  -  Ef 
zm^  21  2um 


^j~*~(B~  -  2B..  E)e'27t  u*f  +  mnu>2u*E 

niQ  2  1  -  0  m 


a(t) 


(5-9) 


H  -  eTtE 


op 


(5-10) 


Equation  (5-9)  yields  the  same  energy  operator  expectation  values  listed  in 
Eq .  (4-2).  It  is  now  straightforward  to  obtain  the  second-quantized  form  of 
physical  quantities  such  as  position,  momentum,  dipole  moment  operator,  etc. 

Let  us  find  the  relations  between  the  time  -  independent  operators  a^(t), 
a(t)  and  time -independent  operators  a^,a.  The  latter  ones  are  simply  defined 
for  a  simple  harmonic  oscillator  as 


a  - 


J 


mo“o 

2X 


x  + 


m0w0 


1 


at  /  ^0 

J  2K 


-LP 

mowo 


Substitution  of  the  above  definitions  into  Eq .  (5-3)  yields 


a(t)  =  ^1(t)a  +  #2 ( t ) a^  +  ^3(t) 


a^(t)  -  0*(t)a  +  0*(t)af  +  0*(t) 


where  the  coefficients  are  given  by 


^(t)  - 


af)u( t)  r  2B  (t)1 

'  1  +  — 4— J 


P2(t) 


aou( t) 

~7T~ 


o 


2B, (t) 


•1  + 


illil 

2 

nr  -* 


07  =  -u(t)B,(t) 


= 


rmoV 


(5-11) 


(5-12) 


(5-13) 


0  "  l  X 


We  now  find  an  explicit  expression  for  the  transition  probability  between 
two  states  |m(t)>  and  |n(0)>,  where  |n(0)>  is  a  simple  harmonic  oscillator 
state  given  by 


|n(0)>  -  |0(0)> 


/n< 


(5-14) 
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Here  |0(0)>  is  the  well-known  ground  state  of  the  simple  harmonic  oscillator. 
Using  Eqs .  (5-6),  (5-12)  and  (5-14),  we  get 

<m(  t)  |  n(0)>  =  y=y  <0(t)  |  (a(t)m(ar)n|0(0)> 

-  7^7  <0(t)  I  (^la  +  V*  +  ^3]m(a+)n|°(0)>  ,  (5-15) 

where  |0(t)>  is  the  time -dependent  ground  state  of  the  damped  driven  harmonic 

oscillator,  given  by  Eq .  (5-8).  <m(t)|n(0)>  represents  the  probability 

amplitude  of  finding  the  system  in  a  state  |m(t)>  at  time  t,  given  that  the 

system  was  in  a  state  |n(0^>  at  time  zero.  To  calculate  the  right  side  of  Eq . 

(5-15),  we  need  to  rearrange  the  operators  in  normal  order,  which  means  that 

all  the  a^'s  appear  to  the  left  of  the  a's.  The  normal -ordered  arrangement  is 

(23) 

carried  out  by  introducing  a  normal -ordering  operator  N  such  that 

[Pxa  +  02a*  +  03lm(af)m  -  n[[o*  +  ,  (5-16) 

where  a  and  a *  are  complex  c-numbers.  With  this  method,  we  obtain 


(5-17) 


kl'k3 

Evidently  a  |0(0)>  vanishes  for  all  k^'s  except  that  .  We  have 

used  above  the  notation  0^  -  /^(t),  0^  -  0  (t)  and  0^  =-  which  will  be 

used  hereafter.  |  j  means  binomial  distribution. 


means  binomial  distribution. 


Substitution  of  Eq .  (5-17)  into  Eq .  (5-15)  gives 


r  i  m-k, 

[m,nl  1 


“>!■«»» -3=fer  5  5  grgg^2^1"2 


k  k 

1  2 

n-k  +k, 

x  <0(t)  |  (a"*")  1  g0(0)> 


(5-18) 


where  [m,n]  -  min(m,n)  implies  the  smaller  of  the  two.  Using  Eqs.  (5-8)  and 
( 5  - 14)  ,  we  find 


n-k.+k9 

<0(t)  I  (a)  g0(0)> 


r  “o  in  rmow  e  $  +  iik  .  *. 

L  n-k1+k2J  [  ,Wl  {  -  ij  eXpl'B3> 


T.*1  "W['K  +  “fK  +  B2X] 


(5-19) 


To  integrate  the  right  side  of  Eq.  (5-19),  we  introduce  a  generating  function 
for  the  Hermite  polynomial  He(y)  such  that 


He(y)  -  Pj  exp{-s2  +  2sy}]s_0  . 


(5-20) 


Then  we  have 


sj^dx  exp[.{B*+-f}x2  +  B*x]  Hn.ki+^(a0x) 


r  J  -I 

-  “7  I(S)  s-0 

LH  c  J 


(i  -  n-k.+k-) 


(5-21) 


where 


e-s2  r  -  r  a2  +  2B*  ?  2Qqs  +  B* 

-  dy  exp  -  - 5 -  y  +  - a -  y 

an  J -m  L  QA 


(5-22) 


which  is  easily  integrated.  Substituting  the  result  of  integration  into  Eq. 
(5-21)  and  performing  some  calculations,  we  obtain  the  following  result: 


*  2 

(B0r 


n-k,+k„ 


r 

*  Ltt 


2  * 
0 +  2Bi 


1  \  rn'ki+k2i 

k  t 

3* 

r2Bi  ■  «oi 

J  2  l  k  J 

k3-0 

(k3/2)! 

k  *  “0 

:  even) 

(5- 

By  Eqs .  (5-18),  (5-19)  and  (5-23),  we  obtain  the  probability  amplitude  as 


*  2 
(B,r 


<m(t)|n(0)>  -  -J—  [-2— « M*1  exp j" - B* ( t )  +  -  2  1 

y2nm!n!  ^  +  2B*J  L  ™1  ?  L  3  2(c2+2B*)J 


[m,nj  m‘kl  n'kl+k2  ,  ,  ,  kl  k2  m-krk2 

r»  l  rkn  fn  i  rn-ki+k2i  <V!  V  v  h 

<k  2-kl)/2 


i  \  5  ywfcjrv) 

_  A  1,  _ A  _ A  X  C  L  J 


k^— 0  k2-0  k^-0 


3  J  (k3/2)! 


r2Bi  -  v 

k3/2 

r  2aoB2  i 

*  2 

*  2 

L2B,  +  oft 

2B.  +  a.-* 

n-k^+k2-k3 


(k^ :  even) 


(5-24) 


The  transition  probability  between  the  DDHO  state  |m(t)>  and  the  simple 
harmonic  oscillator  state  |n(0)>  is  given  by 

Pm  n(t)  -  <m(t) | n(0)> | 2  .  (5-25) 

6  .  Results  and  Discussion 

We  have  modeled  a  molecule  adsorbed  on  a  dielectric  solid  surface  as  a 
damped  harmonic  oscillator  driven  by  a  time -dependent  electric  field, 
consisting  of  a  primary  and  a  secondary  field.  The  Hamiltonian  of  this 
modeled  system  is  not  identical  to  the  mechanical  energy  operator,  showing  a 
general  characteristic  of  the  Kalidola-Kanai- type  Hamiltonian.  The  propagator 
of  Eq.  (3-4)  has  the  same  structure  as  that  of  Gerry ^2^  when  there  exists  no 
driving  force,  and  is  similar  to  that  of  Dodonov  et  al.  '  The  wavefunction 
of  Eq.  (3-7)  has  no  definite  parity,  which  can  be  easily  seen  from  the  matrix 
elements  of  x  and  x2 , 


<m(t)  |x|n( t)>  -  /n+I  u(t)  <5^  +  Jn  u*(t) 

<m(t) |x2 |n(t)>  -  J (n+1) (n+2)  u2(t)  6 m  n+2  + 
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m,n- 


1 


+  E(t)S 


m,n 


7n(n- 1)  u* 


(6-1) 


2 


m,n-2 
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(2n+l) |u(t) I2  +  E2(t) 


m ,  n 


+  2/O^iy  u(  t)E(t)  5m  n+1  +  2 yH  u*(t)E(t)  n.1 


(6-2) 


If  i>  (x,c)  and  \f>  (x,c)  have  definite  parity  and  both  have  the  same  (even  or 
m  n 

odd)  parity,  <m( t) | x | n( t )>  must  be  always  zero,  while  if  they  have  opposite 

2 

parities,  <m(t)|x  |n(t)>  must  vanish.  Equations  (6-1)  and  (6-2),  however,  do 

not  satisfy  this  rule.  This  property  of  the  wavefunctions  does  not  result 

from  the  decay  characteristic  of  the  system,  but  from  the  interaction  between 

the  system  and  the  driving  force,  because  the  wavefunction  of  a  damped  (not 

( 21 ) 

driven)  harmonic  oscillator  has  a  definite  parity.  Due  to  the  absence  of 

definite  parity  of  the  wavefunction,  the  diagonal  element  of  the  dipole  moment 

takes  the  form  given  in  Eq .  (4-5),  in  contrast  to  a  simple  harmonic 

oscillator.  The  real  part  of  the  off-diagonal  elements  of  the  dipole  moment 

(25) 

is  the  same  as  that  of  the  classical  case. 

2 

Because  of  the  time -dependence  (x,t)|  ,  the  expectation  values  of  E 

n  op 

(or  ft)  and  other  physical  quantities  do  not  stay  constant  in  time.  The 
diagonal  parts  of  E  and  ft  are  approximated  as, 


<E  > 
op  n ,  n 


(n  +  |)Xw  e'7t  + 


sliouin 


2m0^wm  ‘  w' 2  +  1 


<H> 


n.n 


-  e7t  <E  > 


op  n ,  n 


(6-3) 


(6-4) 


The  first  term  in  Eq.  (6-3)  shows  the  decay  of  the  quantum  state,  and  the 
second  term  represents  the  energy  absorption  from  the  external  field;  these 
processes  are  related  to  the  lineshape  in  optical  phenomena.  When  there 


2 


exists  no  driving  force  (g (u',S)  -  0),  Eqs.  (6-3)  and  (6-4)  reduce  to  those  of 

(26) 

a  damped  harmonic  oscillator.  Hence,  we  can  say  that  for  a  damped 

harmonic  oscillator  the  expectation  values  of  H  stays  constant  in  time  and 

those  of  E  do  not,  while  for  a  damped  driven  oscillator  the  expectation  of 

both  operators  vary  in  time.  Of  course,  for  the  simple  harmonic  case,  the  two 

operators  are  identical  and  their  expectation  values  remain  stationary. 

The  method  adopted  in  Sec.  5  to  calculate  P  (t)  is  different  from  that 

m,n 

(20  27) 

of  others,  '  and  Eq .  (5-24)  is  of  a  new  form.  In  practice,  however,  it 
is  a  formbidable  task  to  get  the  transition  probability,  from  Eq . 

(5-24).  Even  for  the  simplest  two-level  case  such  that  |m(t)>  =  | 1 ( ; ) >  and 
|n(0)>  -  | 0 ( 0 ) > ,  the  calculation  is  quite  tedious.  Hence,  we  simply  write  the 
result  for  this  case: 


Pl,0  “  l<Kt)|0(0)>| 


o  22..  yt/2 

2qqE  (t)  e' 


1  -  e7t  +  2e"27t  +  (1+  e'")  ft(t) 


r2 
3 


1  +  e 


■7t 


-3 


x  exp 


q2E(t)2  (l  +  e7tC3(t)|i 


1  +  e 


*7t 


(6-5) 


where  we  have  used  the  approximation  D  *  exp(7t/2) ,  since  in  a  real  system 

a>  ™  u>q  and  7  «  w . 

Here  we  note  that  the  transitions  should  occur  on  a  time  scale  shorter 

(23)  -1 

than  the  system  damping  time  7  ,  i.e.,  -yt  «  1.  Thus  Eq.  (6-5)  can  be 

approximated  as 


P1  Q(t)  -  Kt’.t)  exp[ -I(w' ,  t)  ] 


(6-6) 
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I(w*  .t) 


?  9 

9  ^  9  9 

cos  (cj  '  t  -  )  +  (~)  sin  (w't- 

V  1  u) 

O 

'w' 

2mo 

[jw2  -  U'2  +  72/4p  +  (7W')2 

•  (6-7) 


The  transition  probability  depends  strongly  on  the  external  driving  force  and 

2 

the  energy  transfer  mechanism.  Around  the  resonance  frequencies  u'  =  u  + 

2  2 

7^/4  -  w>  ,  Eq .  (6-9)  becomes 


Vw'>  - 


2  2 

«0g  ("'  J) 


2m: 


2  ,2  2..\2  ,  ,  2 
w  -  w  +  7  /4  ^  +  (7W  ) 


(6-8) 


where  this  form  is  related  to  the  lineshape  in  optical  phenomena.  Application 
of  the  results  in  this  work  to  optical  phenomena  on  a  dielectric  surface,  to 

/10\  / O  Q \ 

finding  the  coherent  states^  and  to  the  time -dependent  invariant'1 


problem  will  be  left  for  further  study. 
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